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ABSTRACT 

This  paper  is  concerned  with  the  responses  of  inelastic  rectang¬ 
ular  plates  to  impulsive  loadings.  The  effects  of  strain  hardening  and 
strain  rate  sensitivity  of  the  material  are  taken  into  consideration  in  this 
analysis.  A  variational  principle  in  dynamics  of  inelastic  bodies  subject 
to  finite  deformation  is  used  to  determine  the  deformation  process  of  the 
plate.  A  sandwich  plate  idealization  is  employed.  The  accuracy  of  the 
numerical  solution  is  evaluated  by  comparing  it  with  the  existing  analyti¬ 
cal  and  experimental  results.  The  results  indicate  that  this  method  is 
adequate  for  determining  the  dynamic  behavior  of  inelastic  rectangular 
plates  at  relatively  large  deformations. 
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INTRODUCTION 


The  ^general  problem  of  inelastic  deforrhations  in  structures  subjected 
to  dynamic  loading  has  received  wide  interest  in  recent  years.  There  is 
need  to  correctly  estimate  the  damage  inflicted  as  a  result  of  blast  loadings, 
earthquajkes,.  etc.  ,  on  structures.  These  impulsive  forces  cause  damage 
to  the  structures  in  the  form  of  large  permanent  deformations  w'.ich  usually 
involve  the  effects  of  both  geometric  and  material  constitutive  non-linearities. 

In  most  of  the  previous  studies,  the  concepts  of  limit  analysis  and 
rigid-plastic  idealization  have  been  applied  to  obtain  closed  form  solutions 
for  permanent  plastic  deformations  of  beams,  plates,  etc.  [1-8].  However, 
the  rigid-plastic  idealization  is  applicable  only  to  problems  in  which  the 
elastic  energy  of  a  structure  is  negligible  compared  with  the  plastic  energy 
absorbed  by  the  structure.  Prager  and  Hopkins  [1]  were  among  the  first 
to  investigate  the  load  carrying  capacities  of  a  circular  plate  made  of  per¬ 
fectly  plastic  material  which  obeys  the  yield  condition  of  Tresca  and  the 
associated  flow  rule.  They  assume  small  displacements  and  neglect  mem¬ 
brane  forces.  The  impulse  is  in  the  form  of  a  uniformly  distributed  load 
which  is  applied  suddenly  and  removed  after  a  certain  time  interval.  How¬ 
ever,  later  studies  [2]  have  indicated  that  membrane  forces  do  dominate 
over  the  bending  moments  during  the  final  stages  of  deformation.  In  their 
experimental  work  Griffith  and  Vanzent  [3]  have  shown  that  the  load 
carrying  capacities  of  circular  membranes  are  increased  for  large  inti  ii- 
sity  loads  of  short  duration. 

Wang  and  Hopkins  [4]  have  also  neglected  memi)ran<  Inr*  i-s  in  tlwi  r 
analysis  of  circular  plates.  Pezyna's  |5]  investigation  rc*veal.s  that  the 
pulse  shape  has  very  little  effect  on  the  final  deformation  pattern  of  a 
plate.  Boyd  [6]  has  used  a  numerical  technique  to  solve  the  resulting 
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differential  equation  of  dynamic  deformation  of  a  circular  membrane.  Boyd’s 

analysis  was  based  on  the  deformation  theory  of  plasticity,  lie  therefore, 

neglects  the  effects  of  any  possible  unloading  of  the  material  during  the 

deformation  process.  The  permanent  central  deflections  obtained  by  Boyd 

in 

[6]  and  Witmer  et.  al  [7]  are  found  to  be^^ close  agreement.  This  seems  to 
indicate  that  the  rigid  plastic  assumption  is  adequate  when  membrane  stresses 
are  included  in  the  analytical  model. 

In  all  the  earlier  investigations  reviewed  above  only  the  effects  of 
either  bending  moments  or  membrane  forces  were  considered.  Recently 
Jones  [8]  has  examined  the  combined  effects  of  membrane  forces  and 
bending  moments  on  the  behavior  of  circular  plates  loaded  dynamically  and 
for  the  case  of  large  displacements.  However,  a  simplified  yield  condition 
involving  no-interaction  between  force  and  moment  was  employed  by  Jones. 
Florence  [2]  has  investigated  experimentally  the  behavior  of  simply  sup¬ 
ported  circular  plrtes  subjected  to  uniformly  distributed  impulses.  Ilis 
results  give  ample  evidence  of  the  fact  that  membrane  forces  play  a  very 
important  role  in  large  deformation  of  circular  plates. 

As  the  rigid-plastic  idealization  is  not  appropriate  for  strain  hardening 
materials,  a  more  realistic  method  of  analysis  is  required.  In  this  spirit 
Witmer  et  al  [  9'11  ]  have  extensively  studied  the  axisymmetric  responses 
of  beams,  rings,  circular  plates  and  shells  subjected  to  time  dependent 
loads.  Their  finite  difference  method  of  solving  the  dynamic  equations 
account  for  elastic,  perfecHy  plastic,  strain  hardening  and  strain  rate 
behavior  of  the  material.  However,  because  of  the  small  lime  interval  that 
must  be  used,  considerable  computing  time  is  required  oven  on  a  fast 
digital  computer  in  solving  a  problem.  Cox  and  Morland  [12]  have  dis¬ 
cussed  the  effect  of  a  uniformly  distributed  rectangular  pressure  pulse 
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on  a  simply  supported  square  plate.  Their  analysis  neglected  elastic  strain, 
work  hardening  and  strain  rate  effects.  To  simplify  the  mathematical  anal¬ 
ysis,  a  modified  Tresca's  yield  criterion  was  used  by  them.  Recently 
Jones  et  al  [13]  experimentally  investigated  the  effect  of  uniformly  distri¬ 
buted  velocities  on  fully  clamped  rectangular  plates. 

It  is  well  known  that  variational  principles  have  been  widely  used  for 
solving  static  problems  of  infinitesimal  and  finite  plasticity  [14-17].  A  par¬ 
allel  minimum  principle  in  dynamic  plasticity  has  been  developed  by 
Tamuzh  [18]  for  rigid-plastic  bodies  involving  infinitesimal  deformation. 

It  was  noted  by  Lee  and  Ni  [19]  that  the  principle  is  based  on  a  concept  of 
Gibbs  and  Appell  [20],  in  classical  mechanics,  in  employing  finite  variations 
in  accelerations  in  formulating  a  minimum  principle.  Lee  and  Ni  [19|  have 
advanced  the  concept  to  establish  an  absolute  minimum  prim-iple  in  dynaiim  s 
of  elastic-plastic  continua  subject  to  finite  deformations.  The  minimum 
principle  has  the  advantage  of  circumventing  some  of  the  difficulties  in 
treating  the  loading  and  unloading  material  response.  This  minimum 
principle  is  employed  to  solve  the  titled  problem. 

The  present  study  is  concerned  with  the  analysis  of  the  dynamic  response 
of  a  rectangular  plate  subjected  to  impulsive  loading.  The  impulse  is 
applied  instantaneously  at  zero  time  and  then  decreases  linearly  with  Lime. 
The  response  depends  on  the  non- conservative,  strain  rale  S(;nsitivj!y  and 
strain  hardening  properties  of  the  material.  A  numerical  procedure  liased 
on  the  minimum  principle  is  developed  to  determine  the  maximum  and  the 
permanent  displacements  of  a  rectangular  plate.  The  minimum  principle, 
which  is  derived  from  the  basic  equation  of  motion,  is  written  in  terms 
of  Piola-Kirchhoff  s  stresses,  accelerations  and  Lagrangian  strains. 
Kantorovich’s  method  [21]  is  used  to  determine  the  time  depc.ndeni 
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deformation  profile  of  the  plate.  In  conserving  computing  time,  a  sand¬ 
wich  plate  idealization  is  em  iloyed  in  the  numerical  procedure, 

KINEMATICS 


Because  the  deformed  configuration  due  to  large  deflections  differs 
considerably  from  the  initial  configuration,  the  necessity  of  specifying 
whether  the  stresses  and  strains  are  measured  with  respect  to  the  original 
configuration  or  with  respect  to  the  deformed  configuration  arises.  If  the 
original  state  of  the  material  is  homogeneous  and  stress  free,  it  is  con¬ 
venient  to  express  the  dynamic  equations  in  terms  of  the  original  state.  The 

strain  tensor,  E  referred  to  the  initial  configuration,  is  called  Lagrangian 
u 

strain  tensor  and  is  often  referred  to  as  the  strain  tensor  in  Lagrangian 
coordinates.  Analogously,  the  strain  tensor,  referred  to  the  deformed 
configuration  is  known  as  Almansi's  strain  tensor.  In  the  following,  we 
shall  choose  the  Lagrangian  descriptions  of  stress  and  strain. 

Consider  both  the  original  and  the  deformed  configurations  of  a  body. 

Let  two  systems  of  coordinates  X,  (I  =  1,  2,  3)  and  x. (i  =  1,  2,  3)  be  chosen  to 
describe  the  initial  and  the  deformed  configurations  respectively.  If  we  use 
the  same  rectangular  cartesian  coordinate  system  to  describe  both  the 
original  and  the  deformed  configurations,  then  the  Lagrangian  strain  tensor 
E^jmay  be  defined  in  terms  of  the  displacements,  =  x^-X^  ,  such  that 


^IJ"  2  ’^L,! 


(1) 


Here,  a  partial  differentiation  of  a  variable  with  respect  to  designated 

as  (  ),  K. 

Now  consider  a  rectangular  plate  of  length  2a,  breadth  2b,  and  thick- 
ness  h.  Let  the  components  of  the  displacement  vector  U  of  any  particle 
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originally  located  at  (X,  Y,  Z)  along  the  three  mutually  perpendicular  direc¬ 
tions,  be  u,  V,  and  w  in  the  direction  normal  to  the  plate  respectively.  It 
is  assumed  that  the  plate  is,  thin,  the  magnitudes  of  the  deflections  are  of 
the  same  order  as  the  plate  thickness,  h.  Also,  it  is  assumed  that  the  dis 

placements  vary  linearly  through  the  thickness,  and  is  independent  of  the 
z 

distance^from  the  middle  surface.  Hence,  the  displacement  components  at 
any  point  may  be  expressed  in  terms  of  the  displacements  and  their  deriva 
tives  of  the  middle  surface  as 


u  =  u  -  zw 


(2) 


V  =  V  -  Zw 


in  which  u  and  v  are  the  displacements  along  the  directions  x  and  y 
respectively  of  any  point  on  the  middle  surface  of  the  plate.  Then,  the 


second  time  derivatives  of  the  strain  components  are 

Evv=  u  +(u  u  -iu  +v  +v  v+w  w-fw  ) 

XX  ,X  >X>X  tX  ,X  >x,x  ,x,x  ,x 

••  II.  ^  _  II  ^  —  I!  »  ^ 

E=v  t(u  -fu  u  +v  +v  v  •tw  w  +w  ) 

YY  ,y  ,y  ,y  ,y  ,y  ,y  .y  ,y  ,y  .y 


1 


{3j 


E  =  -:r(U  -iv  -fu.  u  +2u  U  +11  u 
XY  2  '  ,  y  ,x  ,x  ,y  ,x  ,y  .  x  .y 


•  • 


+  2  V  ..  v  +  V  ..  V  +■  v  v  +  2  w  w  _  t  w  w  w  ^ 


,x  ,y  ,  X  ,y  ,  x  ,y 


X  ,  y  .  X  ,y 


Here  a  dot  above  any  term  indicates  partial  differentiation  with  respei  t 
to  time.  The  strain-accelerations  given  by  Eq.  (3)  will  be  employed  in 
the  following  minimum  principle. 


MINIMUM  PRINCIPLE 


Consider  a  body  of  a  continuum  occupying  in  its  natural  stale  a  region 


Hero.  Iho  suom  ripl  s  \,^,z  ’  ro  referred  to  X,  V,/,  •  »..r<l,  .  ?>. 
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and  bounded  by  a  piecewise  smooth  surface  A.  The  body  is  subjected 
to  Ume  dependent  body  force  (per  unit  mass)  and  Lagrangian  surface 
traction  (per  unit  initial  area)  over  that  part  of  the  initial  surface  area 
A^.  At  time  t,  let  be  the  displacement  vector  of  a,  particle  of  the 

body  which  has  an  initial  position  of  (Xj^jin  a  recwngular  Cartesian  co¬ 
ordinate  system.  The  displacements  are  prescribed  over  a  part  of  the 
boundary  surface,  A^.  The  deformation  of  the  body  may  be  described  in 
terms  of  the  Lagrangian  strain  tensor,  defined  by  equation  (1). 

Furthermore  the  Lagrangian  strains  may  be  expressed  as  the  sum  of 
two  parts:  elastic  strain,  E^^  ,  and  plastic  strain  It  is  postulated 

that  the  constitutive  relationships,  in  terms  of  Piola-Kirchhoff  stress  ten¬ 
sor,  S.,,  ,  may  be  strain  velocity  dependent  but  are  not  influenced  by 
strain  accelerations.  In  other  words,  it  is  assumed  that 


^KL  -  ®K1,  EilN-  KiW  e  ’  ''•* 

where  S  is  the  temperature  and  is  velocity  rate  of  plastic 

straining.  The  Piola-Kirchhoff  stresses  satisfy  the  boundary  conditions 


^KL  ^  ^M,  L^  ^T 


(5) 


where  N,,  is  the  outward  anit  normal  to  A  and  5. is  the  Kronecker  sym- 
K  Mi, 

bol. 


D^U 


M 


It  has  been  shown  [19]  that  the  true  acceleration  field,  U  ^ - j 

^  Dr 

of  the  body,  which  has  known  or  predetermined  displacemerl  and  velocity 
fields  at  time  t,  is  distinguished  from  all  kinematically  admissible  on<?8 
by  having  the  minimum  value  of  the  following  functional 
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=  I 


s„,  d  V  +  i  r  p  ij- ,  dv 

y  KL  KL  0  Zjy'^o  M  o 

o 


(o) 


-  I  T, ,  U- .  dA  ~  r  p  F- .  u"  d  V 

J  .  M  M  J,,  M  M  o 

*T 

where  p^  is  the  initial  mass  density.  The  minimum  principle  is  valid  for 
continuous  as  well  as  sectionally  discontinuous  acceleration  fields. 
Ordinarily,  it  is  sufficient  to  use  the  first  variation  with  respect  to  the 
acceleration,  6  J  =  0,  to  establish  governing  equations  or  to  solve  a 

3CC 

problem  by  a  direct  method  of  variational  calculus. 


CONSTITUTIVE  EQUATIONS 


The  stress- strain  relationships  derived  here  are  for  those  materials 
which  are  isotropic,  inelastic  and  work  hardening.  Furthermore  il  is 
assumed  that  the  material  behaves  isothernnally.  It  is  known  that  for  the 
mechanical  behavior  of  the  material  in  fhe  inelastic  range,  the  state  of 
stress  or  strain  may  be  represented  by  a  point  in  a  nine  dimensional  stress 
or  strain  space  respectively.  A  system  of  loading  is  considered  as  a  path 
in  this  stress  space  and  the  corresponding  deformation  history  as  a  path  in 
the  strain  space. 

A  basic  assumption  is  made  that  there  exists  a  scalar  function,  i ailed 
a  yield  function  or  loading  function  F,  which  depends  on  the  states  of  stress 
and  strain  and  the  history  of  loading.  The  equation  F  0  represents  < 
closed  surface  in  the  nine  dimensional  stress  space.  Inelastic  defurni.iiioii 
can  occur  only  when  F  0.  The  condition  F  •;  0  indicates  tluil  no  permanent 
deformation  IS  possible  and  the  condition  F>0  has  no  given  physical  meanings. 
It  is  relatively  easy  to  determine  the  yield  stress  for  a  material  subjeiled 
to  an  axial  tensile  load.  Therefore,  it  is  desirable  to  express  mathematically 
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a  general  stress- strain  relationship  in  terms  of  uniaxial  mechanical  pro¬ 
perties. 

In  a  simple  tensile  test,  work  hardening  means  that  the  stress  is  a 
monotonically  increasing  function  of  increasing  strain.  For  a  general  state 
of  stress,  Drucker's  definition  of  work  hardening  implies  that  energy  cannot 
be  extracted  from  the  material  in  the  process  of  applying  and  removing  a 
system  of  forces  acting  on  a  body.  The  following  conditions  are  a  conse¬ 
quence  of  Drucker's  hypothesis;  a)  The  yield  surface  is  convex  ,  and  b) 
the  plastic  strain  increment  vector  is  normal  to  the  loading  surface. 

Using  these  conditions,  Lee  and  Murphy  [1^  have  obtained  a  relation¬ 
ship  between  Lagrangian  strain  increments  and  Piola-Kirchhoff  s  stress 
increments  as 

E'lj  =  ° Hr, C,  ®kl-  <’> 


and 

E",,  =0  for  F<0  or  dFsO  (8) 

l«l 

2 

where  F  =  J2  "  K  , 

Here  G  is  a  scalor  proportionality  function  which  depends  on  the  cur¬ 
rent  stress,  plastic  strain  and  the  loading  history.  For  a  strain-rate 

2 

insentive  material,  K  is  a  constant  which  is  equal  to  the  maximum  value  of 
J_  that  has  occurred  until  the  current  deformation  state.  J  is  the  second 

z  z 

invariant  of  the  stress  deviation  tensor  S^.  with 


•^2  ""  2  ®  'kL 

^KL  ^  ^KL.  ■  1  ^MM  ^KL 


(9) 


The  expressions  dF>0.  dF=  0  and  dF-^'O  imply  respe<  li  ve| /, 
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neutral  loading  and  unloading  of  the  material  in  relationship  to  the  yield 
surface,  characterized  by  the  loading  function.  For  a  given  state  of  stress 
on  the  yield  surface,  either  loading,  unloading  or  neutral  loading  takes 
place  according  to  whether  stress  increment  vector  is  directed  outward, 
inward  or  along  a  tangent  to  the  loading  surface.  Available  experimental 
results  indicate  that  the  Mises'  yield  function  or  loading  function  leads  to  a 
good  prediction  of  the  initial  yielding  of  an  isotropic  material.  Furthermore, 
in  the  absence  of  Banschinger' s  effect  and  in  a  case  of  relatively  small  strain 
the  Mises'  function  give  also  a  good  prediction  of  the  subsequent  yielding. 
Usually,  the  Mises'  function  has  been  expressed  in  terms  of  the  true  stress. 
For  a  case  of  small  strain,  the  difference  between  the  true  and  Piola- 
Kirchoff  stress  tensors  is  relatively  small.  For  the  thin  plate  problem, 
it  is  expected  the  foregoing  conditions  prevail.  Employing  the  concept  of 
isotropic  hardening  and  assuming  that  the  scalar  function  G  is  a  function 
of  ^2  G  may  be  determined  from  the  results  of  a  simple  tension  test. 

It  may  be  shown  that  [14] 

-r) 

and  (10) 

G  =  0  for  ^2  ^  ® 

where  E^.  and  E  are  the  tangent  and  elastic  moduli  of  the  material 
respectively.  Furthermore  the  general  state  of  stress  is  related  to  the 
simple  tensile  stress,  S,  by  the  relation 

J2  =  S^/3  (11) 

In  simple  tension,  yielding  occurs  when  S  reaches  the  yield  stress 

In  the  case  of  a  thin  rectangular  plate,  the  rate  of  the  second  invariant 
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of  the  stress  deviation  tensor  may  be  written  in  terms  of  the  Pit  la-  Kirchhoff 
stresses  as 

[^XX  (^^XX"  ^YY  )‘^  ^YY  (^^YY  '  ®XX  ) 


+  6  Sjj-Y  j 

The  elastic  strain  rate  given  by  the  Hooke's  law, 

ir'  -  i±^  S  -  S  6  (13) 

^  KL  E  KL  E  ^MM  °  KL 

where  v  is  Poisson's  ratio.  Employing  Eqs.  (7)  to  (13)  and  the  assump¬ 
tion  of  plane  stress,  the  stress- strain  relationship  for  the  thin  plate  may 
be  written  as 

^XX  "  ^11  ^XX  ^12  ^XY  ^13  ^YY 


*  *  A  C? 

^XY  “  ^21  ^XX"^  ^22  ^XY  ^23®YY 
Eyy  =  ^31  ^XX  ^32  ^XY  Ss  ®YY 


where  for  J2  >0  ^z~  ^ 


Sl  = 

-^+  ^  ( 
E^  9  ' 

.^^XX" 

V 

Syy) 

2Sxy  G 

l^^XX" 

^yy) 

^12  = 

3 

X  r 

V^^XX 

■^YYy 

C 

n 

-^+G 

9 

^21  ^ 

S2 

^22  = 

J±Y  X 
E 

2SxY 

G 

2SxY*^ 

(  2Sy  Y  ■ 

■^xx) 

^23  "" 

3 
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Si  ^  Ss 


Sz  ”  S3 

^  _  1  .  ^  (^Sy  "  Sx) 

S3  -  ^  +  G  - ^ - 


for  ^2  ~ 


^11 ""  E’  ^12 ^13^  ~ir  * 


^21  "  ^12’  ^22  "  e'^’  ^23  ~ 

^31  "  ^13’  ^32  ^  ^23’  ^33  ~  4 


In  the  use  of  the  minimum  principle,  it  is  convenient  to  express  stress 
rates  in  terms  of  strain  rates.  Thus,  the  inverted  form  of  equation  (  14) 
may  be  given  as 


Sx  “  ^11  Sx  ■'■  ^12  Sy  ■*■  ^13  Sy 


Sy  ”  ^21  Sx  ^22  ^XY  ■*■  ^23  ^YY 


where 


Syy= 

^31  ^X 

^32 

Sy 

Si  = 

(^22  ^33  ■ 

^22 

;  ^32 

^12  "" 

■(^12^33" 

S3 

^32) 

S3  = 

(^12  ^23  ■ 

Ci3 

^22) 

^21  ^ 

S2 

^22  " 

(^11  ^33  ■ 

^13 

Si) 
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^23  =  -  ( 


*^11 ‘'23'  ‘^13‘'21 


^31  ""  ^13 


^32  “  ^33 


""  (^11^22’  ^12*^21  )^° 


^11  [^22^33~  ^23^32/  ~  ^12  {^21^33^  ^23^3l) 

(19) 

■*■  ^13  (^21^32"  ^31^2) 

In  dynamic  plasticity,  the  mechanical  stress- strain  properties  of  some 
materials  are  signifi.cantly  affected  by  the  rate  of  straining.  The  immediate 
effect  of  straining  rate  is  that  the  yield  stress  is  effectively  increased.  In 
order  to  take  into  account  the  effects  of  the  strain  rate  on  the  yield  stress, 
the  following  strain  rate  formula  which  relates  the  yield  stress  at  any 
instant  of  time  to  the  static  yield  stress  will  be  employed  [11]. 


s'^  =  s  (1+1 

o  o  \  ' 


A  e/t^  t  ll/p 


|i/p) 


where  S‘^  is  the  yield  stress  at  any  instant  of  time  and  S  is  the  static 
yield  stress.  P  and  D  are  material  constants  determined  by  experi¬ 
ments  [10].  e  is  the  "effective"  strain  rate,  and  can  be  expressed  in 


terms  of  the  generalized  strain  rates  as 


•  /  4  \  /'  2  *2  •  •  *2' 

®  (  I  )  \^XX  ^XY  ^XX  ^YY  ^Xy) 


,1/2 
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DERIVATION  OF  FUNCTIONAL 

For  the  case  of  a  thin  flat  plate  subjected  to  an  uniformly  distributed 
transverse  impulsive  load,  by  neglecting  body  forces,  equation  (6)  reduces 
to  the  following  form 

J=[  -^(u  +-  +w^)dV^-J  p(t)waA^ 

Vq  Aj, 

(22) 

(^XX  ^*XX  ^  ^XY  ^XY  ^YY  ^YY  j  ^ 
o 

■'/here  p(t)  represents  the  uniformly  distributed  transverse  load  applied  to 

the  plate.  It  is  seen  that  the  minimum  principle  is  expressed  entirely  in 

terms  of  the  current  state  of  stress,  strain-acceleration  and  acceleration 

not 

of  the  body.  The  current  state  of  stress  is^only  a  function  of  position  but 
also  depends  on  the  loading  history.  This  makes  the  analytical  integration 
of  equation  (22)  very  difficult  to  overcome.  This  difficulty,  an  approxi¬ 
mate  solution,  the  Kanfcororich' s  method  may  be  employed  or  the  functions 
w  (x,  y,  t),  u(x,  y,  t)  and  v  (x,  y,  t)  may  be  represented  by 


m  =  1,  2,  3  . . .  M 
n  =  1,  2,  3,. .  .  N 

Here  X^(x)  Y^(y),  cPj^(x)>  and  Y^(y)  are  chosen  coordinate  functions,  each 
of  which  satisfies  prescribed  geometric  boundary  condlions  and  any  other 


- 
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special  characteristic  of  the  solution  such  as  a  condition  of  symmetry. 
The  coefficients  d  (t),  e  (t)  and  f  (t)  are  unknown  functions  of  time 
t  to  be  determined  by  the  minimization  of  the  functional  J  with  respect 
to  general,  to  obtain  an  exact  solution, 

an  infinite  number  of  coordinate  functions  are  needed.  However,  in 
practice  only  a  finite  number  of  these  coordinate  functions  can  be  con¬ 
sidered.  This  number  depends  not  only  on  the  magnitude  of  the  computa¬ 
tion  involved,  but  on  the  desired  accuracy  of  the  solution. 

IMPULSIVE  LOADING  OF  AN  INELASTIC 
FLAT  PLATE  WITH  EDGE  CONSTRAINTS 

For  verification  the  foregoing  approach  is  applied  to  the  analysis  of 
the  dynamic  response  of  an  inelastic  rectangular  plate  subjected  to  im¬ 
pulsive  loading.  If  I  is  the  applied  impulse  per  unit  area  of  the  plate, 
it  can  be  approximated  by  a  simple  triangular  time  history  starting  at 
pressure  p^  and  reaching  zero  at  a  time  t  such  that  the  area  under 
pressure-time  curve  is  equal  to  this  applied  impulse.  In  other  words 


p  (x,  y,  t)  =  0 

t  <0 

p  {x,y,  t)  = 

tsT 

(24) 

p  (x,  y,  t)  =  0 

t>T 

where  p(x,  y,  t)  is  the  pressure  at  any  time  t.  The  tangential  accelera¬ 
tion  components  u  and  v  of  a  point  away  from  the  middle  surface  are 
expressed  in  terms  of  the  middle  surface  components  u  and  v  by 
differentiating  equation  (2)  twice  with  respect  to  time  they  are 
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Furthermore  the  spatial  coordinate  functions  in  equations  (23)  are 
assumed  to  be  continuous  and  differentiable.  By  differentiating  equation 
(2)  tv/ice  with  respect  to  time  and  once  with  respect  to  space  coordinates 
one  obtains 


u  =  u 

z  w 

,X  ,x 

,  XX 

u  =  U  - 

xw 

>y  >y 

»xy 

V  =  V 

z  w 

,  X  ,  X 

,xy 

V  =  V  - 

z  w 

,y  ,y 

,yy 

(26) 


Using  equations  (3),  (25)  and  (26),  we  can  express  J  given  by  equation 
(22)  in  terms  of  displacements  and  their  time  derivatives  as 


«  I 

J 

0  , 


p  (t)  w  dA.j, 


+  [  is(u  u  +u 

J  y  L  xx'  ,x  ,x  ,x 


X 


o 

2 


+  V  +  V  V 

+  w  w  +  w 

,  X  ,  X 

,  X  >  X  ,  X 

+2S  (u  +  V 

• 

•  • 

h  u  u  +  2u 

xy'  ,y  ,  X 

>  X  ,y  ,  X 

•  «  * 

•  ••  •* 

+  u  u  +  2v 

V  +  V  V  +v 

,X  ,y 

,  X  ,y 

,  X  ,  y  ,  X  ,  y 

+ 

2Ay  w 

+  w  w 

+  w  w  )  +  S 

,  X  ,  y 

,  X  ,  y 

»  X  ,  y'  ^ 

-  2  .  - 

2  -  - 

u  +  u 

u  +  V 

+  V  V  +  w 

y.y  >y  >y  >y  >y  .y 


+  w  ^  )  ^  dV 
>  y  o 


(27) 
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The  volume  integral  J  may  be  expressed  as  a  function  of  acceleration 
coefficients  by  substituting  equations  (23)  and  (25)  into  the  equation  (27). 

The  variation  of  J  with  respect  to  u.  implies  the  variation  of  J  with 

•  «  •* 

respect  to  the  acceleration  coefficients  d,  e  and  f.  The  variation  of  J 
with  respect  to  each  coefficient  results  in  an  expression  which  when 
equated  to  zero,  results  in  the  following  three  sets  of  quasi-linear  algebraic 
equations. 


f  n  ;  P  cp  Y  -  z  d  X  Y  )  (-  zX.  Y.) 

J  L  Z.  ^  o  I  mn^m  '  n  mn  m,  x  n'  '  i,  x  j 


V  m  n 
o 


+  (f  X  tji  -  z  d  X  Y  )  (-  z  X.  Y.  ) 
mn  m’^n  mn  m  n,  y  i  Y 


+  d  X  Y  X.  Y.  ;  ]  dV 
mn  m  n  1  j  ,•  J  o 


X.Y.  dA 
1  J  o 


m  =  1,2,3,... 
n  **  1,  2,  3, ... 
i  r  1,  2,  3, ...  m 
j  =  2,  3, . n 


o 


(28) 


r  '"y  ;  P  i(e  CD  Y  -zd  X  Y  )  cp.  Y.j  1  dV 
J  L  Z  Z.  *^0  I  mn^m  n  mn  m,  x  n  jJJ  c 


V  m  n 
o 


+  ('  I. dV  =  0 
J  ij2  o 


(29) 


Myyp  i(f  X  tj;  -  z  d  X  Y  )X.  ill.  fldV 
o  LLL  o  L  mn  m  n  mn  m  n,  y  i  j  J  J  o 


V  m  n 
o 


t  \  I..,  dV  =  0 

•i  y  O 


(30) 
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where 


1.. .  =  fs  {  -z  X.  Y.  (1  +  u  , 

Ijl  L  XX  i.  i,xx  j  '  ,x 

+  X.  (  -  2  y.  V  +  w  Y.)} 

i,x  j,y  ,x  ,x  yj 

+  S  I  -  2  X.  Y.  (2  +  u  +  V  ) 

xy  I  i,x  j,y  ,x  ,y 

-  2  X.  Y.u  -  2  X.  Y.  V  +  w  X.  Y. 
i>xx  j  ,y  1  j,yy  ,x  ,y  i,x  j 

+  w  X.  Y.  1  +  S  (  -  z  X.  Y. 

,x  1  j,y  J  yy  L  1  j,yy 

(1  +  v  )-zY.  u  X.  +  Y.  w  X. 

,y  J.y  ,Y  i,x  J,  y  ,y  iJJ  (31) 

1.. .  =  i  S  Icp.  Y.  (1+  u  )  [ 

ij2  XX  1^1,  X  j  '  ,  x'  J 

S  |cc.  Y.  (1  +  u  )  +  cp.  Y.  u  , 
xy  I  - 1  J,  y  ,x'  '^1,  X  J  ,  y  J 

+  S  (u  to.  Y.  )  I  (32) 

yy  .y  j»y  J 

1.. -  =  rs  (v  X.  ii(.) 
aj3  L  XX  ,x  i,x^j' 

+  S  X.  (1  +  V  )  +  V  X.  i};.  \ 

xy  1;  i,x  '  ,y'  X  1  '^jjy  ; 

+  S  -j  X.  ij).  (1  +  V  )  )  I  (33) 

yy  1  j»y  .y  J  J 

The  solutions  of  the  quasi-linear  algebraic  equations  are  obtained  by 
employing  a  finite  difference  scheme  together  with  numerical  step-by-step 
method.  The  plate  is  idealized  into  a  number  of  load  carrying  layers 
separated  by  a  core  material.  The  stress  variation  through  each  layer 
is  assumed  to  be  negligible.  The  oniy  function  of  the  core  is  to  maintain 
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an  equal  separation  of  the  layers.  There  is  an  even  number  of  layers  and 
they  are  of  equal  thickness.  The  spacing  between  the  discrete  layers  may  be 
obtained-hy  assuming  both  the  idealized  model  and  the  actual  plate  exhibit 
equal  elastic  bending  stiffnesses.  Thus 

^  ^  "zr"  ^  ^  ^ 

b  =  k  =  2  (34) 

where  k  denotes  the  number  of  layers  and  2b  is  the  distance  between  the 
two  extreme  layers.  The  volume  integrals  in  Eqs.  (28)- (30)  are  reduced  to 
surface  integrals  by  integrating  them  first  with  respect  to  the  z  coordinate. 
The  integration  along  the  z  coordinate  may  be  performed  by  multiplying 
the  integrand  by  a  factor  h/k  and  summing  the  integrand  over  the  number  of 
layers  considered. 

The  surface  integrations  are  carried  out  as  follows.  Let  (M-1)  and 
(N-1)  be  the  number  of  segments  along  the  x  and  y  axes  respectively.  M 
and  N  then  represent  the  number  of  stations  along  the  two  coordinate  axes. 
The  coordinates  (x,  y)  of  any  station  (p,  q)  are  given  by 


P  —  1>  2j  ...  hi 

q  =  1,  2, ...  N 


(35) 


where  a  and  b  represent  the  semi-lateral  or  lateral  dimensions  of  the 
plate,  depending  upon  whether  the  origin  of  the  coordinate  system  is  at  the 
center  or  at  the  corner  of  the  rectangular  plate.  If  f(x,  y)  denotes  the  ' 
integrand  in  a  surface  integral,  the  integration  is  performed  by  using  the 
following  expression 


THEMIS-UKD-71-5 


-  cu  » 


r  ah 

j  f(p,  q)  =  _|2__  f  f(i,  1)  +  F  (1,  NJ  +  f  (M,  1) 

Arp 

qsN~l 

+  f(M,N)+2^  f  (M,q)  1 

q=2 

p=:M-l  p=M-l  q=N~l 

+  2^  (Hp,  1)+  f(p,  N)"*  +  4^  ^  f(p,  q)j  (36) 

p=2  p=2  q=2 

By  the  solutions  of  the  quasi-linear  algebraic  equations,  the  a«  <  eleralions 
and  hence  the  displacement  rates  are  determined.  The  strain  and  stress 
increments  can  be  calculated.  The  new  displacements,  stresses  and 
strains  are  then  computed.  A  time-displacement  history  is  maintain*  d 
throughout  the  motion  of  the  plate  until  it  comes  to  rest  and  then  the  per 
manent  displacements  of  the  plate  are  recorded. 

In  this  analysis  the  detailed  computations  are  carried  out  for  the  lol- 
lowing  two  cases: 

a)  All  edges  of  the  plate  are  simply  supported. 

b)  All  edges  of  the  plate  are  fixed. 

The  boundary  conditions  for  the  first  case  are  satisfied  by  the  following 
chosen  coordinate  functions  (the  center  of  the  plate  is  taken  as  the  origin 
of  the  coordinate  system). 


V  = 
n 


sin 


nny 

b 


■  m  =  1,2,3... 
n  =  1,2,3... 
m  =  1,2,3.. 


n  =  1 , 2,  3  . . . 


(37) 
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It  is  to  he  noted  that  the  functions  X  ,  Y  ,  X  ,  Y  ,  cp  and  f 

m  n  m,  X  n,  y  n 

have  certain  properties.  Substitution  of  Eqs.  (37)  into  Eqs.  (28),  (29)  and 
(30j  yields  the  coefficients 


d  , 


e  . 
U 


P(t) 


4ab 


(2n- l)2m- l)ab 


(.l)"’'*'".  t  L.  iV 

K  ' 

o 


ik=k  I  2  Z  7  !  ^1 

2^’’  1,2,.  2tt  ,,2  i.aiJ 


Ik^l 


■  V/jj  “''o 
o 


he 


ab 
o  4 


h 


] 


and 


- I,.  dV 
‘  V  ‘>>3  ° 


1“  1,  *.•>  nn 
j  a 


(38) 


Ik-k 

where  the  symbol  indicates  the  summation  with  respect  to  each  of 

Ik=l 

the  load  carrying  layers  located  at  z. 

For  the  case  of  the  clamped  plate  with  one  of  the  plate  corners  taken  as 
the  origin  of  the  coordinate  system,  the  edge  conditions  are  satisfied  by  the 


coordinate  functions 
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The  function  X  and  Y  are  not  orthogonal.  However,  the  functions 
m  n  “ 

x’  ^n  y’  '^n  orthogonal.  Now  Eqs.  (28),  (29)  and  (30)  by 

utilizing  Eqs.  (39)  may  be  reduced  to  the  following  expressions 


Ik=k 


ik[z  1 


f  2  • ; 
z  d. 


L. 

m 


in 


TT 

a 


m 


Ik=l 


Y  .  + 
nj 


2  •; 
z  d 


mj 


X. 


im 


+  h  d  X 


mn  mi 


'nil 


-  4p  (t)  ab 


+  *  I..  dV  =  0 

o  ^ 


\ 


p  h  c.  a  Y  .  +  I  . .  dV  =  0 
o  in  nj  .!  y  o 

n  o 


z 


and 


p_f  _  ,  b  X,_„  +J  I., 


L  ■  o  mj 
m 


im 


V  ^j3 


i  =  m 

j  =  J .  .  .  .  ,  n 


(40) 


where 


p2a 

X  .  =  X  X.  dx 

mi  J  mi 

o 


2b 

'  .  =  \  Y  Y.  dy 

nj  J  „  n  J  ^ 


(41) 


The  acceleration  coefficients  d..,  e. .,  and  f. .  are  obtained  by  solving 

ij  ij  iJ  1 

the  three  sets  of  (mxn)  simultaneous  equations  (40). 


RESULTS  AND  DISCUSSION 

Extensive  numerical  results  have  been  obtained  for  a  thin  flat  square 
plate  (a  =  12in.  ,  h=  1/8  and  3/8  in.)  made  of  61S-T6  aluminum  alloy.  The 
stress  strain  relationship  may  be  approximated  by  the  bilinear 
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representation  with  E  =  10  x  10^  psi,  =  40,  000  psi,  and  E^=  10  x  lo"^  psi. 

-  6  2  4 

The  initial  mass  density  is  =  253  x  10  lb/ sec  /in  . 

In  the  analysis  of  inelastic  simply  supported  beams,  Balmer  and 
Witmer  [  11  ]  showed  that  the  maximum  differences  between  the  responses 
for  2  layered  and  4  layered  beams  was  about  8%.  In  the  case  of  a  clamped 
beam,  the  differences  were  slightly  larger.  However,  in  this  analysis 
practically  identical  inelastic  responses  were  obtained  for  2  and  4  layered 
plates.  Moreover,  without  loss  of  occuracy,  it  is  assumed  that  the  in-plane 
displacements  are  infinitesimal. 

For  the  selection  of  the  mesh  size,  a  number  of  trial  runs  indicated 
that  12  and  10  stations  adequately  represent  the  half  span  of  a  simply  sup¬ 
ported  and  a  clamped  square  plate  respectively.'  The  value  of  At  has  been 
chosen  on  the  basis  of  numerical  experiments  and  is  taken  to  be  equal  to 
25  microseconds.  The  input  impulse  considered  here  for  1/8"  and  3/8" 
thick  plates  are  0.  083  and  0.  190  psi- sec.  respectively.  Experimentally 
these  respective  impulses  were  obtained  due  to  a  blast  caused  by  detonation 
of  8,  35  and  36  lbs  charges  of  pentolite  located  centrally  7.  7  and  9.  7  feet 
from  the  face  of  the  square  plate.  The  observed  peak  pressures  were  283 
and  615  psi.  These  impulses  are  approximated  by  simply  triangular  time 
histories  starting  at  the  peak  pressures  and  reaching  zero  at  a  time  such 
that  the  associated  impulses  are  equal  to  the  values  cited  earlier. 

In  Figures  1  and  2,  the  respective  response  time  history  for  simply 
supported  and  clamped  square  plates  are  presented.  In  Figure  1 
the  responses  of  a  simply  supported  circular  plate  by  two  other  analyses 
[6,*^  and  the  permanent  experimental  central  deflection  |  7  ]  arc  also 
shown.  Because  Boyd's  analysis  [6  J  did  not  include  linear  unloading, 

Trail  runs  also  indicated  that  a  nine  term  series  in  equation  ;  23)  was 
adequate  to  represent  the  acceleration  of  motion. 
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the  deformation  process  was  assumed  to  stop  when  maximum  strains 
occurred  in  the  material.  In  Figvire  1,  the  horizontal  line  drawn  from  the 
maximum  displacement  indicates  this  effect.  The  maximum  displacements 
reported  by  Boyd  and  the  present  analysis  are  in  close  agreement  but  differ 
from  those  obtained  by  Witmer  [7].  The  permanent  displacements  predicted 
by  this  analysis  are  in  much  closer  agreement  with  the  experimental  results 
[7]  than  those  of  the  other  two  analyses.  Both  Boyd  and  Witmer  assumed 
perfectly  plastic  material  in  their  investigations. 

From  Figures  1  and  2,  it  is  observed  that  the  plate  reaches  an  equi¬ 
librium  position  soon  after  it  passes  through  its  maximum  displacement. 

By  that  time  the  elastic  energy  has  been  completely  recovered.  The  mean 
vibration  amplitude  following  the  maximum  displacement  are  entirely  in¬ 
elastic.  The  plate  now  oscillates  about  its  new  equilibrium  point  until  it 
finally  comes  to  rest.  The  responses  of  the  center  of  a  simply  support*  d 
and  a  built-in  plate  are  almost  identical. 

Figures  3  and  4,  present  the  deflection  pattern  for  the  two  different 
edge  conditions  considered  here.  Initially  it  is  seen  that  the  deformation 
profile  is  in  the  form  of  a  trapezoid  with  rounded  corners.  This  confirms  ihe 
initial  trapizoidal  profile  which  is  generally  assumed  for  displacements 
in  the  limit  analysis  of  plates.  Any  further  motion  of  the  plate  changes  the 
deformation  profile  to  a  parabolic  shape.  As  the  plate  approaches  its 
permanent  deformed  shape,  the  deformation  profile  seems  to  follow  tlx-  rtome 
shape. 

Deflection  contours  for  different  time  intervals  arc  showi  10  Fa'  in 
5  through  7.  Contours  marked  (a)  to  (c)  arc  for  the  simply  supported  plat* 
and  those  marked  (d)  to  (f)  are  for  the  built-in  plate.  These  contours  arc 
far  apart  in  the  beginning  of  the  motion  but  they  become  closer  as  the  motion 
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proceeds.  The  circuldr  shape  of  these  contours  indicates  the  similar  be¬ 
havior  of  a  square  plate  of  side  2a  and  a  circular  plate  of  diameter  2a. 

In  Figures  8  and  'K  the  membrane  froces  for  the  simply  supported 
plate  along  the  aeini-span  for  t  250  microseconds  and  lor  t  -  750 
microseconds  are  plotted.  It  is  seen  that  the  membrane  forces  are  rtcgiigiWe 
at  the  beginning  of  the  motion  and  become  rather  effective  as  the  motion 
proctu-ds.  That  is,  initially  the  bending  modes  dominate  over  the  membrane 
forces  and  these  roles  are  gradually  rt  versi  d  as  the  plate  i  ome  to  rest. 

in  table  1,  the  permanent  deflections  nl  a  V'  x  5~l/i.’  plate  for  variDUs 
applied  impulses  an  « ompare.-d  with  those  determined  experimentally  by 
Jones  The  results  of  the  present  anilysis  an  lower  than  thus*  by 

Jones,  riie  same  difference  is  observed  iu  the  plots  of  the  perniarn  nt  dt  - 
flections  in  Figures  10  and  12.  Ihis  difference  may  be  attribi  ted  I  *  '.aruiu-. 
factors.  Muwever,  it  is  believed  that  the  major  lactor  may  be  that  the 
experimental  set-up  might  not  precisely  provide  the  fixed  and  .  lamperl  *  dgt 
boundary  conditions  utilized  i  >  the  presint  analysis.  It  is  also  lu  be  i>otc«J 
that  both  the  theoretical  and  experimental  results  indicate  that  th«  piutt.!* 
of  the  permanent  dcTleciion  ol  a  square  plate  su'». acted  to  i^niform  impulH'Vi 
pressure  has  the  shape  of  a  truncated  pyramid. 

CONCLUDING  REMARKS 

A  numerical  procedure  based  un  an  ibHuIote  numnunu  priui  ipi»  «. 
been  ueveloped  for  investigating  the  dv>  amis  respunses  .j  f  ».o, 

gular  plates  at  finite  deformation.  l  iie  minimum  print  q.li  u.i  i  st  tu 
concep  of  limte  variation  in  act  cleratmus  cvIim  h  ir<  un  •  t  i>»  .  .  »  rtei' 

diffic nines  in  treating  the  stress- *^1  rain  r*l  UMinHhip5>i  tiuuui.it  r*  u  >>  . . 

other  principles  involving  virtual  displ.o  tme  nts  ur  '.eh*'  t’lt  n. 
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The*  present  approach  is  a  general  one  since  it  takes  into  account  the 
effect  of  strain  hardening  behavior  and  the  strain  rale  sensitivity  of  the 
material,  WTiich  has  also  the  advantage  of  employing  the  Kantorovich's 
rru  thod  jn  that  a  t  umparati  vel  /  large  time  iinrement  may  be  used  without 
(a<  etiug  lumerical  >iisfahii>ties  in  tiu  infegralHUi  process,  llio  accuracy 
uf  it<e  pre,-ent  nuno  rstal  resolts  l>as  Lt  i  n  ver»fi<  il  by  «.  onipari sion  with  the 
^■^^i^l5^g  aijuiytua],  nurnenial  end  oy;)!  rmo'iiioi  r«  s-idiv  r.bldij.efl  lor  plates 
iinder  inipuisivt  i'  . dines. 


1 


THEHIS-UNI)-71-'> 


-  n  - 


i 


REFERENCES 

1.  Prager,  W.  ,  and  Hopkins,  H.  G.  ,  "On  the  Dynamics  of  Plastic  Circular 
Plates,"  ZAMP  (Journal  of  Applied  Mathematics  and  Phvsics^,  Vol.  5, 
No.  4,  1954,  pp.  317-329. 

2.  Florence,  A.  L.  ,  "Circular  Plate  Under  a  Uniformly  Distributed  Im¬ 
pulse,  "j£terrmtJ££^_Joiirnal_of^olidj_^jrKl_j2l£U£Ujre£,  Vol.  2,  1966, 
pp.  37-47. 

3.  Griffith,  J.  ,  and  Vanzeni,  H.  ,  "Large  Di formation  of  Circular  Mem¬ 
branes  Under  Static  and  Dynamic  L;>adi'ig,‘'  First  International  Con¬ 
gress  on  E.xperimental  Mechanics,  New  York,  196l,  pp.  99-109. 

4.  Wang,  A.  J.  ,  and  Hopkins,  M.  G.  ,  "On  the  Plastic  Deformation  of 
Built  in  Circular  Plates  under  Impulsive  Load,"  Journal  of  Mechanics 
and  Physics  of  Solids,  V  >!.  3,  No.  1,  l'^54,  pp.  22-37. 

5.  Pezyna,  P.  ,  "Dynamic  Load  Carrynng  Capacity  of  a  Circular  Plate," 
Arch.  Mech.  Stos.  ,  V'ol.  10,  No.  5.  19"S,  pp.  o3S-647. 

h.  Boyd,  D.  E.  ,  "Dynamic  Deformation  ot  Circ  uiar  Membranes,  " 

Journal  of  the  Engineering  Mechanics  Divisiu n ,  Proceedings  of  the 
American  Society  of  Civil  Engineering,  Vol.  No,  EM3,  June  Idbo, 

pp.  1-lb, 

7.  Witmer,  E.  A.  ,  Balmer,  11.  A.,  Leech,  J.  W.  ,  and  Pian,  T.  H.  II.  , 

"Large  Dynamic  Deformations  of  Beams,  Rings,  Plates  and  .Shells," 
American  Institute  of  Aeronautics  and  Astronautics  Journal,  Vol.  1, 
r96'3,  pp.  !B48-i85b. 

8.  Jones,  N.  ,  "Impulsive  Loading  of  a  Simply  .Supported  Circular  Plastii. 
Plate,  "Journal  of  Applied  Mechanics,  March  1968,  pp.  59-65. 

9.  Leech,  J.  W.  ,  Witrncr,  E.  A.,  and  Pian,  T.ll.  II.,  "Numerical  Cal¬ 
culation  lechniuuc  lor  Large  Elastic-Plastic  Transient  Deformations 
of  rhin  Shells,"  American  Institute  of  Aeronautics  and  Astronautics 
.lournal ,  Vol.  o.  No.  TTi  1968,  pp.  2352-2359. 

10.  Leech,  J.  W.  ,  "Finite-Difference  Calculation  .M  ethod  for  Large 
Elastic-Plastic  Dynamically-Induced  Deformations  of  General  Thin 
Shells,  "  Air  Force  Flight  Dynamics  Lab,  AFFDL- 'I'R -  66- 1  7 1 , 

December  l9bb. 

11.  Balmer,  H.  A.  ,  and  Witiner,  E.  A.  ,  "Theoretical  and  Lxpe rinient  cl  Cor¬ 
relation  of  Largt  Dynamic  and  Permanent  Deformations  of  Impulsively 
Loaded  Simple  Structures,"  Technical  Documentary  Report  No.  t'DL- 

T  DR  -  64-  108 ,  Air  Force  f'lighl  Dynaniics  Laboratory,  Wrighl-Patterson 
Air  Force  Base,  Ohio,  I9bi. 

12.  Co.x,  A.  D.  ,  and  Morland,  1..  W.  ,  "Dynamic  Plasfu  Deformations  of 
.Simply  Supported  Square  Plate's,'  Journed  of  Mech.  Phys.  Solids, 

Vol.  7.  22u,^>4!. 


THKMIb-UND- /i-b 


~  c.p>  - 


13.  Jones,  N.  ,  Uran,  T.  O.  ,  and  Tekin,  S.  A. ,  'The  Dynamic  Plastic 
Behavior  of  Fully  Clamped  Rect*  igular  Plates,"  Int.  J.  Solids 
Structures,  1970,  Vol.  6,  pp.  i4''*«1512. 

14.  Lee,  L.  H.  N.  ,  and  Murphy,  L.  M.  ,  "Inelastic  Axj  symmetric  Buckling 
of  Ring  Stiffened  Cylindrical  Shells  Usider  External  Pressure," 

Technical  Report,  Thernis-UND-'69-2,  University  of  Notre  Dame,  l‘)69. 

15.  Drucker,  D.  C.  ,  "Variational  Principles  in  the  Mathematical  Thetsry 
of  Plasticity,  "  Proc.  Symposia  in  Appl.  Math.,  Vol.  Vlll.  Calculus 
of  Variations  and  Its  Applications.  McGraw-Hill  Book  Co.  , 

pp.  7-22. 

16.  Sanders,  J.  L.  ,  McComb,  H.  G.  ,  and  Schlechte,  F  R.  ,  ’’A  Variational 
Theorem  for  Creep  with  Applications  to  Plates  and  Columns,  NACA, 
Report  1342,  1957. 

17.  Roth,  R.  S.  ,  "Plastic  Buckling  of  Thin  Shallow'  .Spherical  Shells,"' 
Proceedings  of  the  Fourth  U.  S.  National  Congress  of  Applied  Mechanics, 
ASME,  1962,  pp.  1059-1065. 

18.  Tamuzh,  V.  P.  ,  On  A  Minimum  Prin>.i>)!e  in  Dynamics  ol  Rigid- Plastic- 
Bodies,  "  Pk^  Vol.  26,  No.  4,  1967,  pp.  71S.';22. 

19.  Lee,  L.  H.  N.  ,  and  C.M.  Ni.  ,  "A  Minimum  Principle  in  Dynamics  of 
Elastic- Plastic  Continua  at  Finite  Deformation,"  A.  naper  submitted  to 
the  International  Symposium  on  Foundations  of  f'lastn  ity,  Warsaw, 
Poland,  September,  1972. 

20.  Pars,  L.  A.  ,  "A  Treatise  on  Analytical  Dynamu.s,"  Ih-inernan,  London, 
1965,  p.  200. 

21.  Kantorovich,  L.  V.  and  Krylov,  V.  I.  .  "Approximate -Methods  'jf  Higher 
Analysis,"  Interscience  Publisher,  New  York,  1958  p.  240. 

22.  Sureshwara,  B.  ,  "Impulsive  Loading  of  Rectangular  Plates  with 
Finite  Plastic  Deformations,"  dissertation  submitted  to  the  Graduate- 
School,  University  of  Notre  Dame.  Notre  Dame,  Indiana,  in  August, 

1971,  in  partial  fulfillment  of  the  requirements  for  the  degree  of 
Doctor  of  Philosophy. 


Ni0ll031d30 


DISTANCE  FROM  CLAMPED  END  (IN.) 

Middle  Span  Deflections  of  a  Clamped  Squai'e  Plate 


Deflection  Contours  of  a  Square  Plate 


DISTANCE  FROM  THE  CENTER  OF  THE  PLATE  (IN.) 

Fig.  8  Middle  Span  Membrane  Forces  of  a  Simply  Supported  Square  Plate, 


11  Middle  Span  Permanent  Profile  of  Aluminum  6061  T6  Plate 


